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1. Effective potential as a function of the average plaquette 

The study of the QCD phase diagram at non-zero temperature (T) and chemical potential 
{ptq) is one of the most important topics among studies of lattice QCD. In particular, the study of 
the endpoint of the first order phase transition hne in the {T,}Xq) plane is particularly interesting 
both from the experimental and theoretical point of view. The existence of such a critical point 
is suggested by phenomenological studies. The appearance of the critical endpoint in the {T^ptq) 
plane is closely related to hadronic fluctuations in heavy ion collisions and may be experimentally 
examined by an event-by-event analysis of heavy ion collisions. Many trials have been made to 
prove the existence of the critical endpoint by first principle calculation in lattice QCD, however 
no definite conclusion on this issue is obtained so far. The purpose of this study is to clarify the 
existence of the endpoint of the first order phase transition line in the {T^ptq) plane. We propose a 
new method to investigate the nature of transition. 

We evaluate an effective potential as a function of the average plaquette (P) ' , and identify the 
type of transition from the shape of the potential. The partition function can be written as ^ 

irGc(i8,At,) = j&U {Atm{^iq)t e-'^^''^^^ = j /?(P,At,)w(P)e-^«(^-^) dP, (1.1) 

where Sg{P,p) is the gauge action, w{P) is the state density at /i^ = for each P, 

= w{P',p) = J&U d{P'-P) {dctMf'e^P^^-^^, (1.2) 

and R{P,jJLq) is the modification (reweighting) factor for finite ptq, which is defined by 

J^{^,^q)- j^u 5{P'-P){detM{0))f^i {S{P'-P)) 



(i8.M<,=0) 



where (• • •)(j3 ^^=o) means the expectation value at /i^ = 0, detM is the quark determinant, Nf is 
the number of flavors ^ , and A'^^ite = A'^^^ x Nt is the number of sites. We then define the effective 
potential as V {P, p , jJ-q) = — ln(/?we^'^s). If there is a first order phase transition point, where two 
different states coexist, the potential must have two minima at two different values of P. In this 
paper, we discuss whether the potential at /i^ = 0, i.e. ln(we^'^^), which is quadratic function when 
the transition is a crossover, can change to a double-well potential by the reweighting factor for 
finite Hq, as illustrated in Fig. [T] (left). 

2. Taylor expansion in terms oijiq/T and Gaussian distribution 

However, the calculation of the quark determinant is quite expensive and is actually difficult 
except on small lattices. Moreover, the calculation of R{P,ii.q) becomes increasingly more difficult 



'For later discussions, we define the average plaquette P as P= —Sg/ (6j3A'site)- This is the average of the plaquette 
over all elementary squares for the standard gauge action. 

^We restrict ourselves to discuss only the case when the quark matrix does not depend on j8 explicitly, e.g. the 
standard Wilson and staggered quark actions, the p4-improved staggered quark action etc., for simplicity. 



A'f must be replaced in these equations to Nf/4 when we use a staggered type quark action. 
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Crossover Re weighting factor 




e=(A'/4)Im[ln(det Af)] 

Figure 1: (left) Schematic figures of the effective potential and the reweighting factor (right) The histogram 
of the complex phase for ^q/T = 1.0 and 2.0 at /3 = 3.65. 



for large Hq due to the sign problem, i.e. the statistical error becomes exponentially larger as }iq 
increases. We avoid these problems by the following two ideas. One is that we perform a Taylor 
expansion of IndetM(jLi^) in terms of ji^ at /i^ = and calculate the expansion coefficients [|l|], 

^N'sNttD,\^)\ (2.1) 

The Taylor expansion coefficients are rather easy to calculate by using the stochastic noise method. 
Although we must cut off this expansion at an appropriate order in ptq, we can estimate the appli- 
cation range where the approximation is valid for each analysis [Q, ^]. While the application range 
of the Taylor expansion of In should be limited by the critical point because In 3^qc is singular 
at the critical point, there is no such Umit for the application range in the expansion of \nR{P,pLq) 
because the weight factor should always be well-defined. 

The sign problem is avoided by the following idea. We consider a probability distribution func- 
tion as a function of the complex phase of the quark determinant 0, |F| = | detM(jU^)/detM(0)|^f 
and P, i.e. w{P,\F\,Q). If we assume the distribution function in d is well-approximated by a 
Gaussian function, the sign problem in the calculation of \nR{P,}Xq) is completely solved iQ]. 

We define the complex phase by a Taylor expansion, d = N{N^NtY,'^^QlmD2n+i{tJ.q/TY"'^^ . 
Since the partition function is real even at non-zero density, the distribution function has the 
symmetry under the change from 6 to —6. Therefore, the distribution function is written by 
w{d) '-^ exp[— ((320^ + + a^d^ + ■■■)]• Moreover, because D„ is a trace of a matrix which 
has space index [^], e.g. D\ oc ti[M^^dM/d{lJ.q/T)], the central limit theorem suggests that the 
distribution function is well-approximated by a Gaussian function, when the system size is suffi- 
ciently large in comparison to the correlation length between diagonal elements of the matrix. We 
plotted the distribution of the complex phase in Fig. |l| (right) and fitted by a Gaussian function 
(dashed line). It is found that this approximation is quite well, hence we consider the leading term 
of the expansion only, w{P, \F\,6) oc y^a2{P,\F\)/K exp[— a2(^, 1^1)^^]- The coefficient a2{P, \F\) 



In 



detM{iJ,q) 
detM(O) 



17=1 



5"(lndetM) 
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Figure 2: (left) The plaquette histogram and the effective potential at fiq — 0. (right) The reweighting factor 
ln7;(P,Al<j) for Al^/r = 0.5-2.5. 



is given by 1/ (2a2) = (6 ) for each P and The numerator of Eq. (1.3) is then evaluated by 



{Fin,)5{P'-P)) 



(2.2) 



Because l/ai 0(A^site)> the phase factor in R{P,}JLa) decreases exponentially as a function of the 



volume. However, the operator in Eq. ( |2.2[ ) is always real and positive for each configuration in this 
framework, hence the expectation value of R{P, ptq) is always larger than its statistical error, namely 
the contribution \nR{P,}iq) to the effective potential V{P,^,}JLq) is always well-defined. Therefore, 
the sign problem is completely avoided if we can assume the Gaussian distribution of 6. 



3. Numerical results of the effective potential 

We calculate w{P,j5) and R{P, jXq) using data obtained by simulations in [^. The Taylor expan- 
sion coefficients are computed up to 0{}JLq). These are measured at sixteen simulation points from 
j8 = 3.52 to 4.00 for the bare quark mass ma = 0. 1. The corresponding temperature normalized by 
the pseudo-critical temperature is in the range of T /Tc = 0.76 to 1.98. The ratio of pseudo-scalar 
and vector meson masses is nips /my ~ 0.7 at jS =3.65. The lattice sizeA^siteis 16^ x4. The number 
of configurations is 1000 - 4000 for each j8. 

The probability distribution function w(P,j8), i.e. the histogram of P for each j3, and the 
effective potential j3) at /i^ = are given in Fig. ^left). To obtain w{P,P), we grouped the 
configurations by the value of P into blocks and counted the number of configurations in these 
blocks, and the potential V(P,P) is normalized by the minimum value for each temperature. 

The results for \nR{P, ptq) are shown by solid lines in Fig. ^ (right) for ptq/T = 0.5, 1 .0, 1 .5, 2.0 
and 2.5. We find a rapid change in In/? around P ~ 0.83, and the variation becomes larger as piq/T 
increases. The dashed lines in Fig. ^ (right) are the results that we obtained when the effect of the 
complex phase, i.e. exp[— 1/(4(32)], is omitted. These dashed lines correspond to the reweighting 
factor with non-zero isospin chemical potential jXi and zero quark chemical potential ptq [|5|]. The 
variation of \nR in terms of P becomes milder when the effect of the complex phase is omitted. 
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Figure 3: The slope (left) and curvature (right) of lnR{P, jXq). The dashed line is the curvature of — Inw. 

This explains the difference between the phase diagrams of QCD with non-zero quark chemical 
potential and non-zero isospin chemical potential [Q]. 

We discuss the shape of the effective potential at non-zero /i^. The effective potential is ob- 
tained from V(P, j8,/x^) = — lnw(P, j8) - ln/?(P,/i^) substituting the data in Fig. ^ From Eq. ([T^), 
the slope of V{P,^,lJLq) can be controlled by j3, 



V(P,j8,At,) = V(P, A),Ai,) - 6(j8 - A))A^site/'. 



(3.1) 



under a change of j3o ^ jS, however the curvature of the potential does not change by j3. We expect 
that the curvature vanishes at the endpoint of the first order phase transition line by canceling 
cf{\nw)/dP'^ and d^{lnR)/dP^. In order to analyze the sign of d'^V /dP'^{P,}Xq), we fitted the data 
of InP by a quadratic function of P and calculate the first and second derivatives of lnP(P,jU^) at 
each P. The results of the slope and curvature are shown in Fig. |3| for each piq/T . In the region 
around P ~ 0.83, d{\nR) jdP becomes larger as \XqlT increases and \xyR(P^\iq) changes sharply in 
this region. The magnitude of the curvature of InP also becomes larger as \lqlT increases. 

To evaluate d^{Va.w) j dP^{P), we assume w{P,l5) in Fig. ^(left) is a Gaussian function. In this 
case, the curvature at (P) is given by —d^{lnw)/dP^ = 6A^site/ZP> where Xp = 6Nsite{{P — {P))'^) 
is the plaquette susceptibility. The dashed line in Fig. ^ (right) is the result of —d'^{\nw)/dP'^{P). 

It is found from Fig. ^ (right) that the maximum value of d^ {In R)/dP^{P, jig) at P = 0.80 
becomes larger than — (In w)/<iP^ for Hq/T ^ 2.5. This means that the curvature of the effective 
potential vanishes at Hq/T ~ 2.5 and becomes negative for large ptq/T , namely the shape of the 
effective potential which is of quadratic type at /x^ = changes to a double-well type at large 
piq/T . For the quantitative estimation of the endpoint of the first order phase transition, further 
investigation must be needed. However, this argument strongly suggests the existence of the first 
order phase transition line in the {T,jXq) plane. Further details of this analysis are given in iQ]. 

4. Canonical partition function 



Next, we want to apply the effective potential argument to the weight factor as a function of 
the quark number density p. The physical meaning of this potential is clearer than that of P because 
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Figure 4: The chemical potential as a function of the quark number density. 



the weight factor for each quark number corresponds to the canonical partition function !^c, 



(4.1) 



N 



The canonical partition function can be given by an inverse Laplace transformation ^ § 

3 r""/^ 



3fc{T,N) = ^l ^(At//r)^'-^(«'/^+'^^'/^)irGc(r,Ato + /M/) 

271 J-k/3 



(4.2) 



where fM) is an appropriate real constant. Note that f^Qc{T ,}Xq + 2TiiT /3) = ^Gc{T,fJ.q)- Recently, 
this canonical partition function is calculated for Nf = 4 using the Glasgow method [^. However, 
with present day computer resources, the calculation by the Glasgow method is difficult except on 
small lattices. We consider approximations which is valid for large volume and low density in this 
approach, as discussed in the first part of this paper. 



We calculate the grand partition function by the Taylor expansion, Eq. (2.1), 



1 



f<i-mN)Y\^,,M{0)re-'-^ ^ (ei'^^'y^^-Mn] ) (4.3) 



\ detM(O) 



where V = N^. We moreover use a saddle point approximation, which is valid for a large sys- 
tem. We find a saddle point zo in the complex Hq/T plane for each configuration, which satisfies 



-1 



J^Gc(r,0) exp 



0. The canonical partition function is given by 

-'«/2 , / 



V NfNt £ Dnzl - pzo 



1 



V|D"(zo) 



) (4.4) 



for large W. Here, D"{z) = {d'-/dz^) {NiNti::=iDnZ") and D"(zo) = |Z)"(zo)k'"- 

The chemical potential, i.e. the slope of the effective potential, is also evaluated by 



T 



-\d\n^c{T,pV) 



zo exp [V (NfNt D„z"o - pzo)] e 



1 

V|D"(zo)| 



dp 



exp \y [NiNtYrn=x - pzo)] e- 



V|D''(zo)| 



.(4.5) 



(r,/j,=o) 
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This equation is similar to the formula of the reweighting method for finite jXq. The operator in 
the denominator corresponds to a reweighting factor, and the chemical potential is an expectation 
value of the saddle point calculated with this modification factor. 

We analyze the data used in the previous section. The Taylor expansion coefficients up to 
0{jXq) are used. The volume V = 16^ would be sufficiently large, and we assume a Gaussian 
distribution function for the complex phase of the reweighting factor, again. We find a saddle point 
Zo numerically for each configuration, assuming zo exists near the real axis in the low density region 



of the complex IXq/T plane. We use multi-jS reweighting method [10] combining all data obtained 
at 16 points of j3. Configurations are generated with the provability of the Boltzmann weight in 
Monte-Carlo simulations, however the important configurations will change when the weight is 
changed by the reweighting method. For such a case, the multi-j8 reweighting is effective, since the 
important configurations are automatically selected among all configurations generated at multi-j8, 
and also this method is useful for the interpolation between the simulation points. 

We plot the result of jXq/T in Fig. as a function of p /T^ for each temperature. The dot- 
dashed line is the value in the free gas limit. As seen in Fig. ^ (left), the configurations do not 
distribute uniformly in the range of P which is necessary in this analysis, and correct results cannot 
be obtained if the important configurations are missing. At low temperature, the important value 
of P changes very much as p increases, therefore we plotted only the data when the expectation 
value of P is on the peaks of the histograms of P in Fig. |2| (left). The dashed lines are cubic spline 
interpolations of these data. 

It is found from this figure that a qualitative feature of ^q/T changes around T /Tc ~ 0.8, i.e. 
ptq/T increases monotonically as p increases above 0.8, whereas it shows an s-shape below 0.8. 
This means that there is more than one values of p /T^ for one value of Hq/T below T /Tc ~ 0.8. 
This is a signature of a first order phase transition. The critical value of jXq/T is about 2.5, which is 
consistent with the result in the previous section. Although further studies including justifications 
of these approximations used in this analysis are necessary for more qualitative investigation, this 
result also suggests the existence of the first order phase transition line in the {T,}jLq) plane. 
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